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We report on the fabrication of high-Q, fused-quartz microresonators and the parametric genera- 
tion of a frequency comb with 36 GHz line spacing using them. We have characterized the intrinsic 
stability of the comb in both the time and frequency domains to assess its suitability for future 
precision metrology applications. Intensity autocorrelation measurements and line-by-line comb 
control reveal near-transform-limited picosecond pulse trains that are associated with good relative 
phase and amplitude stability of the comb lines. The comb's 36 GHz line spacing can be readily 
photodetected, which enables measurements of its intrinsic and absolute phase fluctuations. 



Femtosecond-laser optical frequency combs have revo- 
lutionized frequency metrology and precision timekeep- 
ing by providing a dense set of absolute reference lines 
spanning more than an octave. These sources exhibit 
sub- femtosecond timing jitter corresponding to an ul- 
tralow phase noise of < 100 /zrad on high harmonics of 
their typically 100's of MHz repetition frequency (line 
spacing) [TJ. This remarkable level of performance has 
enabled measurements of atomic transition frequencies 
at the 17th digit [2] and the generation of ultralow noise 
signals [3 . Beyond their natural application as an optical 
clockwork, frequency combs are used in diverse applica- 
tions including precision measurements of fundamental 
physics [TJ E] , direct spectroscopy and real-time trace gas 
detection [T], molecular fingerprinting [5], astronomical 
spectrograph calibration |6 , and optical arbitrary wave- 
form generation [7]. An even broader range of applica- 
tions may be possible if the bulk and complexity of a 
femtosecond comb could be reduced without significantly 
increasing fluctuations in the comb spectrum. 

A new class of frequency combs has recently emerged 
based on monolithic microresonators, henceforth denoted 
microcombs [5l412j . Here the comb generation relies on 
parametric conversion provided simply by third-order 
nonlinear optical effects and is enabled by recent ad- 
vances in the quality factor Q and the small volume 
of microresonators. These devices require only a single 
continuous-wave laser source, but the usable frequency 
span of the comb depends on low dispersion, making 
material properties critical. Microcombs also present 
a unique platform for creating large line spacings (10's 
of GHz to THz). To date microcombs possessing 100's 
of lines each spaced by 100's of GHz have been created 
with silica microtoroids [5] and silicon-nitride microrings 
[TUl ITT] . Uniformity [5j and control [T3] of silica micro- 
combs have also been demonstrated. Microcombs with 
mode spacings of 10 to 25 GHz have also been realized 
with machined crystalline resonators |9j and aspects of 
their microwave-domain spectral purities have been ana- 
lyzed 0U|. 

Future frequency metrology applications of micro- 
combs will require a line spacing (repetition rate, / rcp ) in 
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FIG. 1. High-Q fused-quartz microresonators. (a) A fused 
quartz rod is rotated in a lathe (not shown) and rubbed 
against a fixed metal form charged with diamond abrasive. 
The ground form is then flame polished with household-grade 
propane and oxygen (not shown), (b) Image showing an ar- 
ray of fused quartz microresonators after grinding and flame 
polishing. The disks of 1.9 mm diameter are separated by 
0.45 mm. (c) Frequency scan of a mode for disk four (indi- 
cated by the arrow) demonstrating an optical Qq — 5.2 x 10 8 . 
(d) Cavity ringdown at K = 1 to determine the optical Q of 
disk four. The dashed lines indicate a decay time of 260 ns, 
corresponding to Qo = 6.2 x 10 s . 

the 10's of GHz (millimeter-scale resonator), comb span 
approaching an octave, and low absolute phase and fre- 
quency fluctuations that can be further reduced by wide- 
band comb control [T5] . Additionally, a threshold power 
for comb generation in the mW range, and the capability 
for integration with chip-based photonic circuits would 
enable novel portable applications. Here we present a 
platform towards attaining these requirements. Our fab- 
rication approach takes advantage of precision mechan- 
ical shaping and polishing possible with fused quartz, 
a versatile amorphous material with low optical losses 
throughout the visible and into the infrared. We created 
resonators that feature high Q and low effective opti- 
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FIG. 2. Microresonator lineshape and corresponding comb 
spectra for laser pump power of 62 mW. (a) Thermally broad- 
ened lineshape with a mostly sawtooth shape. Power levels 
Li, L,2, L3 respectively correspond to the parametric thresh- 
old, the initiation of a 36 GHz comb but with unstable be- 
havior, and the region of stable 36 GHz comb generation, (b) 
Comb spectrum at Li. (c) Comb spectrum at L2. (d) Comb 
spectrum at L3. 



cal mode area, which allowed comb generation for a <10 
mW pump laser. By pumping above threshold, a comb 
with ~ 150 lines spaced by 36 GHz was generated. We 
harnessed full line-by-line control (amplitude and phase) 
of the comb to generate near-transform limited optical 
waveforms and to understand the comb's relative phase 
stability. In the frequency domain, we characterized the 
comb's mode-spacing fluctuations to assess prospects for 
stabilizing it to optical atomic standards. 

We fabricated an array of disk-like microresonators on 
a common substrate by use of a combination of diamond 
grinding and flame polishing. A schematic of our proce- 
dure is shown in Fig. [IJi. Fused quartz rods with « 2 mm 
diameter were rotated in a ball-bearing spindle. To create 
the basic shape of a disk, the glass was rubbed against a 
metal form with triangular cross section. Diamond abra- 
sive between them removed material from both. At a 
diamond abrasive size of 6 /xm, ~30 minutes of grinding 
was required to create the underlying shape shown in Fig. 
[TJd. The resulting surface was far too rough to support 
high Q whispering-gallery modes. By polishing the disks 
with a propane-oxygen flame we achieved an extremely 
smooth surface. Substantial melting of the quartz was 
required to alleviate subsurface damage caused by grind- 
ing. Figure shows an image of a completed disk ar- 
ray. Free-spectral range measurements of several disks 
revealed a diameter uniformity of 0.1 %. 

After processing, the fused quartz rod was placed in 
a temperature-stabilized mount (<0.01 K stability) and 



whispering-gallery modes of the disks were accessed us- 
ing standard fiber tapers [T5TtT7] . The coupling between 
disk modes and the taper can be characterized by the 
parameter K — tq/t & , where t is the intrinsic photon 
lifetime and r e is the photon loss rate due to coupling. 
The frequency linewidth of the resonator-taper system 
is 7 = Tq -1 + t" 1 . Figure fit shows a frequency scan 
over a disk mode using a tunable diode laser at 1560 
nm. A low optical power of ^10 /iW was used to mini- 
mize thermal broadening of the resonance [18] and taper 
coupling was small (K ss 0.03). The linewidth of the 
Lorentzian transmission feature is 0.37 MHz, which cor- 
responds to an optical Qq — ojt of 5.2 x 10 s , where 
lo is the optical carrier frequency |19j . We confirmed 
this Qq measurement by recording the resonator energy 
ringdown time following a rapid frequency scan across 
resonance with K = 1; see Fig. [Tji. A fit to the cavity 
energy decay yields Qq = 6.2 x 10 8 , in good agreement 
with our linewidth measurement. Material absorption 
in the fused-quartz samples we use imposes a limit to 
Qo ~ 2 x 10 10 . We speculate that some contamination 
occurs during the flame-polishing step |20j . 

To achieve the long-term stability in K needed for 
studies of microcomb generation, we place the taper di- 
rectly in contact with the resonator surface. Translating 
the taper along the axis of the quartz rod allows us to 
tune the resonator-taper coupling from under- to over- 
coupling. Moreover, measurements of K at various lo- 
cations along the disk edge provide a qualitative picture 
of the mode's spatial profile. Specifically, the resonator 
mode studied in Fig. [l]c and d has two field nodes in the 
axial direction and is presumably of lowest radial order. 
Numerical calculations based on an ideal disk model [ST] 
suggest the effective area of this mode (A e g) is approxi- 
mately 500 /im 2 . The threshold for parametric oscillation 
threshold is P th = \ i±f £ £ jfi , where n (n 2 ) is the 
refractive index (nonlinear index) of fused quartz and Ai> 
is the free spectral range. For our disk resonators we ex- 
pect P t h = 4 mW for K = 0.2. 

We have generated a microcomb at, and significantly 
above, parametric threshold. Figure [2^-d shows a disk 
resonator transmission lineshape and typical comb op- 
tical spectra. With constant laser power, we adjusted 
and stabilized the intracavity power by tuning the pump 
laser frequency for a fixed K — 0.2. Here the lineshape is 
thermally broadened and assumes a characteristic "saw- 
tooth" shape [18]. Power level Li indicates the observed 
3.7 mW threshold for parametric oscillation, which is in 
excellent agreement with our prediction. Notably, the 
signal and idler modes are spaced by 828 GHz (23 x 36 
GHz); indicating the point at which resonator- induced 
phase mismatch of pump, signal, and idler balances the 
mismatch created by nonlinear effects j9]|22j[23]. Increas- 
ing the intracavity power further resulted in a widened 
comb span up to ±25 nm about 1560 nm and a reduction 
of the comb spacing toward the fundamental 36 GHz of 
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FIG. 3. Time-domain measurements of a three line comb, (a) 
Schematic of our setup to control and analyze the microcomb 
spectrum. The amplitude and phase of light exciting the res- 
onator via tapered fiber is manipulated with a programmable 
filter. An intensity autocorrelator is used to record the fields' 
waveforms, (b) Optical spectrum following the filter and am- 
plifier, (c) Time-averaged autocorrelation signal for A<p — 
(black circles) and A(j> — ir (blue circles), which respectively 
demonstrate constructive and destructive interference of the 
three comb fields. The solid lines show the results of an an- 
alytic model, (d) Visibility of constructive interference as a 
function of A<j>. 



the disk resonator, which is consistent with the simula- 
tions of Chembo et al. [2H |2SJ • However, our observa- 
tions indicate that 36 GHz comb generation is complex 
with competition between different modes of operation. 
Specifically, we present two instances of 36 GHz combs 
with dramatically different properties: (1) At power level 
L 2 (Fig. [2J3) an «36 GHz comb spacing is initiated, but 
a finer analysis (see Fig. |6fj) reveals numerous comb line 
spacing frequencies and (2) in a range about power level 
L3 (Fig. [2ji) a robust comb pattern exists with a sin- 
gle line spacing of 36.0012 GHz. Interestingly, when we 
tuned to power level L2 a deviation in the sawtooth shape 
of the transmission lineshapc was observed. This behav- 
ior depends on resonator-taper coupling and laser polar- 
ization, but not the laser frequency scan rate, which we 
varied from 0.005 to 1 s. Moreover, in this range stable 
time-domain waveforms and low noise frequency-domain 
performance were not evident under the conditions of 
Fig. [2] [26]. About L 3 we observed the resumption of a 
sawtooth lineshape. 

The spectral structure similar to Fig. [2] has been ob- 
served in a variety of parametric microcombs, but little is 
known about the temporal structure of the output as de- 
termined by the phase relationship between the comb el- 
ements. Measurements and control of optical waveforms 
can reveal crucial information about the internal fields 
of the microcomb [27] [25] ■ Using time-domain diagnos- 



tics, we studied the relative phase stability of the first 
three comb lines (pump, signal, and idler) at paramet- 
ric threshold, and the phase stability of many lines far 
above threshold. For a comb composed of only three 
equal-amplitude lines its time waveform has two distinct 
behaviors, which depend primarily on the relative phases 
of the lines: constructive interference, resulting in a rudi- 
mentary pulse at /rep, or partial destructive interference, 
resulting in a beat signal at 2/ rop . A comb with many 
lines exhibits more complicated behavior and offers high 
peak power waveforms. 

Our experimental starting point was a comb spectrum 
similar to that shown in Fig. [2]d, except with spacing 
of 864 GHz [29]. Using a programable amplitude and 
phase filter featuring 10 GHz resolution (Fig. [3^,), we 
selected the central three lines of this comb, equalized 
their amplitudes to better than 7 %, and adjusted their 
phases. The optical spectrum we obtained following am- 
plification with a low-dispersion erbium fiber amplifier is 
shown in Fig. [3]d. It was delivered via a short section of 
optical fiber to an intensity autocorrelator for character- 
ization. By varying the relative phase (A(f>) of the lowest 
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FIG. 4. Generation of picosecond pulses utilizing line-by-line 
comb control, (a) Measured time-domain waveform (black 
line) following amplitude and phase optimization, and sim- 
ulation (red dashed line) based on the measured amplitudes 
of comb lines used and assuming phase alignment of all the 
fields. The inset shows the optical spectrum after optimiza- 
tion, (b) Optical spectrum following broadening in highly 
nonlinear fiber. The comb was amplified to 1.6 W. 
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wavelength line, we were able to explore the transition 
between constructive and deconstructive interference of 
the comb lines. The relative phases between the center 
and the highest wavelength line were set to zero. Fig- 
ure [3J; shows the time-waveform autocorrelation signal 
(P(t)), where r is the differential path delay, for A<fi = 
(black open circles) and A0 = tt (blue open circles). The 
solid lines indicate the results of an analytic model based 
on three phase-coherent fields, which are in excellent 
agreement with the data. To characterize the interfer- 
ence between the three fields we introduce a visibility V 
based on P(r), which is defined as V — I p(^/a^)+p(o) I- 
The filled points in Fig. [3ji are measurements of V as 
a function of A<j>, and the red line shows our analytic 
model. We measured a maximum (minimum) visibility 
of 0.7 (0.007), and the model predicts the visibility is 
0.768 (0.002). This indicates good phase coherence of 
the comb lines. Fitting the model to the data with free 
parameter A<p yields an uncertainty in A(f> of 20 mrad, 
which we associate with an upper limit on the fluctua- 
tions of the relative phase of the comb elements over the 
measurement time of 1700 seconds |30j . 

We also carried out time-domain waveform generation 
with a microcomb far above parametric threshold. The 
amplitude and phase filter was used to adjust the central 
15 lines (0.54 THz) of a 36 GHz comb similar to that of 
Fig. [2ji. Here a manual line-by-line procedure was used 
to optimize the peak waveform intensity. This resulted 
in stable near-transform-limited 2.5 ps optical pulses, as 
demonstrated by the time-averaged autocorrelation sig- 
nal in Fig. [4^,. Hence, the relative phases of the 15 comb 
lines were presumed to be constant across the spectrum 
and were correspondingly stationary in time. Moreover, 
the autocorrelation signal remained unchanged despite 
large perturbations to the microcomb system, such as 
re-locking the pump laser frequency or re-placing the 
tapered fiber in contact with the disk. These observa- 
tions indicate that the internal fields of the resonator and 
the pulse generation mechanism are deterministic and 
repeatable. Significantly, the creation of time-domain 
waveforms, such as shown in Fig. 4a, implies that such 
a compact and simple device could be competitive with 
more conventional mode-locked lasers a source of high 
repetition rate ultrashort optical pulses that would be 
useful for a variety of time- and frequency-domain appli- 
cations |12) . For example, nonlinear broadening in fibers 
could be a route to ultra-broadband spectra for comb 
self-referencing and comb spectroscopy. Using a 101 m 
length of highly nonlinear fiber, we have broadened our 
15 line comb by more than a factor of 10, as shown in 
Fig. [4Jd. This demonstrates sufficient peak power to drive 
nonlinear processes external to the microresonator. 

Frequency-domain techniques offer significant advan- 
tages in characterizing microcomb fluctuations, includ- 
ing access to the power spectral density, wide measure- 
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FIG. 5. Residual phase noise of microcomb line spacing, (a) 
Schematic of our experimental setup. The microcomb spec- 
trum is split (BS) into two paths, which are separately filtered 
at Ai = 1560 nm and A2 = 1553.6 nm. Photodectors PDj.,2 
convert the comb spacing of the independent paths into elec- 
trical signals near 36 GHz. PDi (PD2) generates 2.15 mA 
(0.04 mA) of photocurrent. Double-balanced mixers Mi and 
M2 convert these signals to baseband, and mixer M3 offsets 
one path by 11 MHz. This offset frequency was referenced to 
a hydrogen maser. At mixer M4 the two paths are interfered 
yielding residual noise, (b) Single-sideband phase noise of 
the microcomb 36 GHz tone (black line) and the contribution 
from the 11 MHz offset oscillator (red line), (c) Histogram 
of zero-dead-time counter measurements of the signal exiting 
M4. (d) Allan deviation of the microcomb line spacing (filled 
points). The solid line shows the 1/r scaling exhibited at 
most averaging times. 



ment bandwidth, and extremely high precision. A key 
feature of our system is a comb line spacing of 36 GHz, 
which enables direct photodetection and analysis of the 
resulting microwave signal in the frequency domain. Here 
we present measurements of the microcomb's residual 
and absolute phase fluctuations. Residual (i.e. intrin- 
sic) noise indicates the degree that the frequencies of 
different comb lines are correlated, and absolute noise 
indicates the stability of the lines with respect to a fixed 
frequency reference. Moreover, we have characterized the 
relationship between line spacing fluctuations and pump 
laser intensity noise. This information provides a bench- 
mark for what level of control will be required for future 
applications of stabilized microcombs. 

To understand the intrinsic stability of the microcomb, 
we directly compared the line spacing at two independent 
portions of its spectrum. We created two copies of the mi- 
crocomb spectrum shown in Fig. [2ji, using a fiber beam- 
splitter, and separately bandpass filtered them at 1560 
nm and 1553.6 nm; see Fig. [5^,. Independent high-speed 
photodetectors created electronic signals at 36 GHz. The 
frequencies of these signals were reduced to baseband 
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FIG. 6. Microcomb fluctuations measured with respect to 
fixed references. For all of these experiments the entire mi- 
crocomb spectrum was delivered to PD2. (a) Free-running mi- 
crowave spectrum of the microcomb operating at power level 
L3 obtained via photodetection. (b) Power spectral density 
(PSD) of L^> (black), line-spacing signal amplitude (red), and 
pump laser RIN (blue). The RIN attains its shot noise lim- 
ited value of -151 dBc/Hz beyond 100 kHz. The cyan shaded 
region shows a scaling of the RIN to phase noise using the 
measured conversion factor 7p for our system. The top range 
of this area corresponds to a constant 7p with frequency, and 
the bottom range accounts for the measured frequency re- 
sponse of pump power to line spacing. The green line shows 
the phase noise contribution of the 35.5 GHz oscillator, (c) 
Microwave spectrum at L2 with a 30 kHz resolution band- 
width. Under these conditions the comb exhibits multiple 
line spacings near 36 GHz. (d) (black) and amplitude 
(red) noise of the microcomb at L2. The blue trace shows 
pump laser RIN. 



with a common 35.1 GHz oscillator at mixers Mi and 
Ma, and the path corresponding to filter Ai was further 
shifted by 11 MHz with a low-phase-noise oscillator. A 
final mixer (M4) interfered the two baseband signals to 
reject common-mode fluctuations. Hence, only residual 
phase noise of the microcomb (and noise associated with 
the independent paths) appeared at the output of M4. 
Using a commercial phase noise analyzer we recorded the 
single-sideband phase noise spectral density L^, which is 
shown by the black line in Fig. [5]d. For reference, the red 
line in Fig. [5]d shows the phase noise contribution of the 
11 MHz offset signal. The microcomb's residual phase 
noise is extremely low. At carrier offsets greater than 10 
kHz, is dominated by photodetection noise and would 
be reduced with greater microcomb optical power. Close 
to carrier the residual phase noise is mostly below -100 
dBc/Hz, a value comparable to the absolute phase noise 
of the best optical and microwave oscillators |3J. This 
indicates that with appropriate frequency control such 
a microcomb could be harnessed for portable low-noise 
synthesis applications. 

We have assessed whether the microcomb line spacing 



is the same at 1560 nm and 1553.6 nm by electronically 
counting the 11 MHz signal exiting M4. Here an offset 
from 11 MHz would indicate non-equidistance of the mi- 
crocomb spectrum. Previous work by Del'Haye et. al. 
studied microcomb equidistance by way of an auxiliary 
fiber-laser frequency comb [8] . With access to the 36 GHz 
comb spacing, our measurements directly characterize 
microcomb equidistance for the first time. We acquired 
546 consecutive (zero dead time) one-second long mea- 
surements; Fig. [5]: shows a histogram of the frequency 
difference (A/) between the microcomb and reference os- 
cillator. The data exhibit a Gaussian distribution with 
mean of -0.65 //Hz with a 25 /iHz width. Hence, the mode 
spacing of the comb does not change fractionally by more 
than 3xl0 -17 across 6.4 nm. The Allan deviation <ta, 
shown in Fig. [5ji, demonstrates the intrinsic stability of 
the microcomb line spacing with cta=1 x 10~ 17 after only 
400 s of averaging. 

Many future applications of a microcomb will require 
a stable spectrum with respect to fixed frequency refer- 
ences, such as atomic clocks. Here we show the micro- 
comb line spacing's microwave spectrum, absolute phase 
noise, amplitude noise, and its dependence on pump laser 
power. Importantly, different behavior of these was ob- 
served at power levels L2 and L3 that was not revealed 
in coarse measurements of the optical spectrum (Fig. [5J; 
and d). The experimental setup for these measurements 
was similar to Fig. [5^i, except that the beamsplitter, op- 
tical filters, and mixers Mi 3,4 were removed. The pho- 
tocurrent generated at PD2 was converted to baseband 
with a 35.5 GHz oscillator for analysis. For phase noise 
measurements we used a digital prescaler (divide-by-265) 
following Mi to reject spurious amplitude fluctuations. 

First, we present our observations at L3. Figure [6^, 
shows a single microwave tone 55 dB above a broad asym- 
metric pedestal. The width of this feature is limited by 
the resolution bandwidth of 10 kHz. The amplitude and 
phase fluctuations of this tone are shown by the red and 
black traces in Fig. [6]3, respectively. Notably, the ab- 
solute phase noise is between 20 and 120 dB above the 
residual noise, depending on the carrier offset. For com- 
parison the green line shows the phase noise of our com- 
mercial 35.5 GHz oscillator. The pump laser power, with 
the relative intensity noise (RIN) spectrum shown by the 
blue trace in Fig. [BJa, significantly influences the comb 
line spacing. We measured a line spacing-power depen- 
dence of |tp| = 0.3 MHz/mW with a 3 dB bandwidth of 
^20 Hz. The correlation between pump laser RIN and 
microcomb amplitude and phase noise is evident in Fig. 
[6]d. Scaling the RIN by jp results in the phase noise 
prediction shown by the cyan-shaded area. Pump laser 
power fluctuations explain a significant fraction of the 
line spacing phase noise, and either passive or active re- 
duction of the RIN will be crucial for future experiments. 
In particular, obtaining a f-2f heterodyne beatnote [3"lT 
133) . a key outstanding milestone in microcomb systems, 
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will be challenging in the face of significant comb noise. 

Figure |6p and d shows a similar analysis of the micro- 
comb line spacing's spectrum (c), amplitude (red trace in 
(d)) and phase fluctuations (black trace in (d)) at power 
level L2. Under this operating mode, the comb exhibits 
multiple line spacings in a few MHz range around 36 
GHz. We observe significantly larger (10's of dB) ampli- 
tude and phase noise, even though the pump laser R1N 
is the same as at L 3 . The correspondence between pump 
RIN and amplitude, phase noise remains evident. We 
have not been able to develop a detailed understanding 
of, nor any strategy to mitigate, the poor performance 
at L2. However, these observations emphasize that mea- 
surements beyond those of the optical spectra will be crit- 
ical to understanding the properties of the microcomb in 
all its operating regimes. 

In summary, we have fabricated mm-scale, high-Q op- 
tical microresonators with fused quartz. By using them, 
a frequency comb with 36 GHz line spacing, 50 nm span, 
and <10 mW threshold was created near 1560 nm. We 
studied the microcomb spectrum in both the time and 
frequency domains. Picosecond optical pulses were gen- 
erated by way of full line-by-line control of up to 15 
comb lines. Precise frequency-domain techniques enabled 
a direct test of comb line fractional equidistance at the 
3xl0 -17 level, and a characterization of the microcomb's 
absolute line spacing phase noise. Quantifying the noise 
of microresonator frequency combs is crucial in assess- 
ing whether they can distribute modern microwave and 
optical atomic standards. In the future we will focus gen- 
erating a wider comb span, and stabilizing the comb to 
Rb transitions at 780 nm by way of higher bandwidth 
comb control. 
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hala and Tobias Kippenberg, and we thank Gabe Yeas 
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manuscript. This work is a contribution of NIST and 
is not subject to copyright in the United States. 
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